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Abstract. We propose the formulation of a dihadron fragmentation function in 
terms of parton matrix elements. Under the collinear factorization approximation 
and facilitated by the cut- vertex technique, the two hadron inclusive cross section 
at leading order (LO) in e+ e- annihilation is shown to factorize into a short 
distance parton cross section and the long distance dihadron fragmentation 
function. We also derive the DGLAP evolution equation of this function at leading 
log. The evolution equation for the non-singlet and singlet quark fragmentation 
function and the gluon fragmentation function are solved numerically with the 
initial condition taken from event generators. Modifications to the dihadron 
fragmentation function from higher twist corrections in DIS off nuclei are 
computed. Results are presented for cases of physical interest. 



PACS numbers: 13.66.Bc, 25.75.Gz, 11.15.Bt 
1. Introduction 

One of the most promising signatures for the formation of a quark gluon plasma 
(QGP) in a heavy-ion collision has been that of jet quenching 1 . This phenomenon 
leads to the suppression of high pt particles emanating from such collisions. Such jet 
quenching phenomena have been among the most striking experimental discoveries 
from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory. 
Simultaneously, the modification of jets in cold nuclear matter has been measured via 
the deep-inelastic scattering (DIS) of leptons off nuclei at the HERMES detector at 
the Deutsche Elektronen-Synchrotron (DESY). Such experiments, besides providing 
an essential baseline for jet modifications in hot matter, are interesting tests of nuclear 
enhanced higher twist effects on parton propagation in nuclei. 

In the investigation of jet suppression in heavy-ion collisions, correlations between 
two highpT hadrons in azimuthal angle are used to study the change of jet structure 
While the back-to-back correlations arc suppressed in central Au + Au collisions, 
indicating parton energy loss, the same-side correlations remain approximately the 
same as in p + p and d + Au collisions. An almost identical situation has been 
observed in the modification of fragmentation functions in cold nuclear matter: while 
the single inclusive measurements indicate a large suppression hadrons with large 
momentum fractions (large z), the ratio of the double inclusive measurements to the 
single inclusive measurements show minimal variation with atomic number. Given the 
experimental kinematics, this is considered as an indication of parton hadronization 
outside the medium. However, since the same-side correlation corresponds to two- 
hadron distribution within a single jet, the observed phenomenon is highly nontrivial. 
To answer the question as to why a parton with a reduced energy would give the 
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same two-hadron distribution, one has to take a closer look at the single and double 
hadron fragmentation functions and their modification in a medium. In this article, we 
present results for the evolution of the dihadron fragmentation function in the process 
of e + e~ annihilation, and its modification in a cold nuclear medium. Modifications in 
a hot QGP will be presented in a later effort. 

For reactions at an energy scale much above Aqcd, one can factorize the cross 
section into a short-distance parton cross section which is computable order by order 
as a scries in a s (Q 2 ) and a long-distance phenomenological object (the fragmentation 
function) which contains the non-perturbative information of parton hadronization [3] . 
These functions are process independent. If measured at one energy scale, they may 
be predicted for all other energy scales via the Dokshitzer-Gribov-Lipatov-Altarclli- 
Parisi (DGLAP) evolution equations j^. In this article, we will extend this formalism 
to the double inclusive fragmentation function D^ 1 ' 112 (z\, Z2, Q 2 ) or the dihadron 
fragmentation function, where Zi,z 2 represent the forward momentum fractions of 
the two hadrons. In the absence of experimental data, the dihadron fragmentation 
functions will be sampled at one energy scale from event generators and the prediction 
of the derived DGLAP evolution will be verified against the measurement at another 
energy scale. The dihadron fragmentation function will then be subjected to medium 
modification in a DIS scenario. Comparisons with data from the HERMES experiment 
will be presented. 



2. The double fragmentation function and its evolution 

We begin our analysis with the following semi-inclusive process 

e + + e~ — > 7* — > hi + hi + X 

of e + e~ annihilation. Wc consider two-jet events where both the identified hadrons hi 
and h 2 emanate from the same jet. At leading order in the strong coupling, this occurs 
from the conversion of the virtual photon into a back-to-back quark and anti-quark 
pair which fragment into two jets of hadrons. 

In the limit of very large Q 2 of the reaction, we may invoke the collinear 
approximation. Under this approximation, at leading twist, we can demonstrate the 
factorization of the two-hadron inclusive cross section into a hard total partonic cross 
section o-q and the double inclusive fragmentation function D^ lh2 (zi, Z2) (see Ref. [S] 
for details), 



D h q ^(zi,z 2 ) + D'^ h2 (zi,z 2 )\ . (1) 

In the above equation, the leading order double inclusive fragmentation function of a 
quark is obtained as 
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In the above equation z — z\ + z 2 , p represents the momentum of the fragmenting 
parton, is the sum of the momenta of the two detected hadrons i.e. pi + p 2 , q± is 
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the relative transverse momentum between the two detected hadrons and n is a light- 
like null vector. The sum over S indicates a sum over all possible final hadronic states. 
The above equation may be represented by the diagrams of the cut vertex notation 
|Sj as that in the left panel of Fig. ^ Note that all transverse momentum q± up to a 
scale have been integrated over into the definition of the fragmentation function. 
Hadrons with transverse momenta > n± may not emanate from the fragmentation of 
a single parton. 




Figure 1. The left panel represents the cut- vertex representation of the dihadron 
fragmentation function. The right panel represents a Next-to-Leading order 
correction from quark and gluon single fragmentation. 

In the interest of simplicity we first specialize to the case of the non-singlet 
fragmentation function, D^g 2 (zi,Z2) = D^ l,hi {zi,Z2) — Dq 1 ' 3 {z\,Z2)- With the 
definition of the dihadron fragmentation functions in the operator formalism, which is 
shown to factorize from the hard parton cross section at LO, we may now evaluate its 
DGLAP evolution by computing the double inclusive cross section at next to leading 
order (NLO). This constitutes a rather involved procedure; the reader is referred to 
Ref. [5] for details. The resulting DGLAP evolution consists of two parts. There 
remains the usual evolution due to the radiation of soft gluons and both detected 
hadrons emanating with a small transverse momentum off the same parton. Hadron 
pairs with transverse momentum larger than fj,± arc produced pcrturbatively, by the 
independent single fragmentation of a quark and a gluon which emanate from the 
splitting of a quark. The contribution of this process to the evolution of the dihadron 
fragmentation function is represented by the cut-vertex diagram in the right panel of 
Fig. ^ Incorporating this diagram, we obtain the DGLAP evolution of the dihadron 
fragmentation function as 



High Pt hadron-hadron correlations 



4 



dD h N t(Q* 
dlogQ 2 



2tt 



*qg * ^NS 



(Q 2 ) + P q 



> qg -*D h N i s (Q 2 )D h g *(Q 2 ) + l 



(3) 



where the switch 1 — > 2 is meant solely for the last term. The expressions for the 
splitting functions ( P q ^ qg and P q ^ qg ), single fragmentation functions, as well as the 
convolution notations (* and *) may be obtained from Ref. [Hj- The above equation 
may be solved numerically and results for the cases of z\ = 2z% and 3^2 are presented 
in Fig. We assume a simple ansatz for the initial condition at Q 2 = 2 GcV 2 i.e. 
D t ^g 2 (z 1 ,z 2 ) = D hl (zi)D h2 (z 2 ). We present results for the evolution with log(Q 2 ) at 
intervals of 1, up to log(Q 2 ) = 4.693 i.e. Q 2 = 109 GeV 2 . It may be noted that the 
results demonstrate minimal change as the energy scale is varied over a wide range of 
values. 
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Figure 2. Results of the ratio of the non-singlet quark dihadron 
fragmentation function D q % 2 (z±,Z2,Q 2 ) to the single leading fragmentation 
function Dq 1 (z±, Q 2 ). In the left panel z\ = 2z2, in the right panel z\ = 3z2- 
Results are presented from Q 2 = 2GeV 2 to 109GeV 2 . 



In generalizing from the non-singlet quark fragmentation function to the quark or 
antiquark fragmentation function one simply includes a piece in which the fragmenting 
quark may radiate a gluon and both detected hadrons may emanate from the 
fragmentation of this gluon. However, this new piece couples the evolution of the 
quark fragmentation function with the evolution of the gluon fragmentation function. 
The evolution of the gluon fragmenation function has the usual components where the 
gluon may radiate further gluons prior to undergoing fragmentation, or convert into 
a quark-antiquark pair with the detected hadrons materializing from either parton. 
It also contains the new additional pieces where each of the detected hadrons may 
emanate from each of the radiated partons. We illustrate the various components of 
the gluon fragmentation function by the following equation (sec Ref. |Zj for details), 
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In the above equation, / indicates the flavour of a fermion. In the first term, the 
sum runs from 1 to 6 indicating a sum over both quarks and antiquarks. In the second 
term the sum runs from 1 to 3 indicating only fermions. The case of quarks replaced 
with antiquarks is indicated with the switch 1 to 2 which only applies to this term. The 
remaining terms indicate splitting into two gluons followed by fragmentation function. 
Expressions for the splitting functions may be obtained from Ref. [7]. 

Once the fragmentation functions are measured at a given energy scale, their 
variation with energy scale can be predicted with the aid of the evolution equations 
mentioned above. To date, there remain no measurements of two particle correlations 
in e + e~collisions. However, Monte Carlo event generators such as JETSET have 
enjoyed great success as simulators of such collisions. The quark and gluon dihadron 
fragmentation at a scale /z 2 = 2GcV 2 are extracted from two and three jet events 
in JETSET (see Ref. for details). These are represented by the filled symbols 
(triangles for quarks and squares for gluons) in Fig. In both cases, the momentum 
fraction of the leading hadron is set as z\ = 0.5. Once measured, these fragmentation 
functions are then evolved with the aid of the evolution equations outlined above. 
Results of the evolution with log(Q 2 ) at intervals of 1, up to log(Q 2 ) = 4.693 i.e. 
Q 2 = 109 GeV 2 are presented in the left panel of Fig. [31 as the various lines. The 
fragmentation functions are once again measured at the higher scale from JETSET 
and compared with the calculation (open triangles for quarks and open squares for the 
gluons). The excellent agreement between the evolution equation and the predictions 
from JETSET at the higher scale provide an important test of the derived DGLAP 
equation for the dihadron fragmentation function. 
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Figure 3. Left panel shows results of the evolution of the quark and gluon 
dihadron fragmentation functions. The right panel is a plot of the ratio of 
D hlh2 (z± , Z2, Q 2 ) to the single leading fragmentation function D hl (z\, Q 2 ), for 
the cases shown in the left panel. See text for details 



The right panel of Fig.Elshows the variation of the ratio D(zi, 22, Q 2 )/D(zi, Q 2 ) 
with Q 2 while z\ is held fixed at 0.5. While the quark dihadron fragmentation function 
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displays a slow variation, the gluon fragmentation function shows almost no change 
as the energy scale is varied over two orders of magnitude. This is consistent with the 
results obtained for the evolution of the non-singlet functions. 

3. Medium Modification 

The HERMES detector at DESY measures the semi-inclusive cross section for the 
production of one or two hadrons from electroproduction off nuclei such as Nitrogen 
and Krypton The measured quantity. R^ , is the number of events with at 
least two detected hadrons with the leading hadron assuming a momentum fraction of 
z\ > 0.5, divided by the number of events with at least one hadron with a momentum 
fraction z > 0.5. 

Due to various detector systcmatics, a double ratio is presented: Ra/d = 
i?^ l/l2 /i?p 1 ' 12 , where i?^, 1 ' 12 is the same ratio as discussed above but measured in 
scattering off Deuterium. Assuming that the Deuterium nucleus is two small to induce 
any noticeable medium effects on the fragmentation of hard partons passing through 
it and in the limit of low pair multiplicity per event, we equate the above ratio as 

A/D D^ h *(Q 2 )/D^(Q 2 y 
where D^^^D^ 1 represent the quark fragmentation functions in vacuum, and 
E)h 1 h2^£)h 1 re p rcscn t those modified by the medium. Values of Ra/d f° r Nitrogen 
are plotted in the right panel of Fig. 0] 

The calculation of the medium modification of the fragmentation function is 
carried out in the nuclear enhanced higher twist formalism. This formalism is 
presented in detail in Ref. In principle, one calculates the cross section for the 
process where the recoiling lepton imparts a large momentum transfer {Q 2 ) to one of 
the quarks in the incoming nucleus via single photon exchange. The struck quark then 
undergoes multiple scattering off the soft gluon field of the nucleus prior to exiting the 
medium and fragmenting to hadrons. In this process it absorbs a certain amount 
of transverse momentum and virtuality, which it loses by gluon bremsstrahlung. 
Alternatively, the struck quark may be produced off shell and may radiate gluons 
prior to scattering off the soft gluons. In cither case the radiated gluon (which also 
carries a colour charge) may also encounter further scattering of the gluon field prior 
to exiting the medium and fragmenting into a jet of hadrons. Interferences between 
the various processes mentioned above also need to be considered. 

The eventual calculation of the medium modification involves the computation of 
multiple higher twist diagrams such as the one displayed in the left panel of Fig. 0] 
While these diagrams, which involve quark gluon correlations in nuclei, are suppressed 
by powers of Q 2 , they receive an enhancement ~ As from the fact that the quark and 
gluon may originate from any of the A nuclcons in the nucleus (for further details see 
Refs. PIE]). Calculations for a Nitrogen nucleus are presented and compared to the 
HERMES data in the right panel of Fig. H 

4. Discussions and Conclusions 

This study was motivated by the observation [2] that the same side correlations of 
two high px hadrons in central Au + Au collisions remain approximately unchanged 
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Figure 4. Left panel shows one of the higher twist diagrams involved in the 
calculation of the medium modification of the fragmentation functions. The 
various dashed lines indicate the various cuts: left (L), central (C), right (R) 
and virtual (V). The right panel shows the experimental data from HERMES for 
Ra/d an d f ne solid line is the calculation. Both the data and the calculation are 
preliminary. 



as compared with that in p + p collisions. A similar situation was also observed in DIS 
experiments. Neglecting the differences in production cross section and fragmentation 
functions for different parton species, the integrated yield of the same side correlation 
over a small range of angles should be the ratio of the dihadron to the single 
hadron fragmentation functions, D l ^ lh2 (zi,Z2,Q 2 )/D , l ! l 1 (zi,Q 2 ). Thus, we evaluated 
the influence of DGLAP evolution on the ratio D% lfl2 (zi, z 2 , Q 2 )/D^ (zi, Q 2 ). Our 
numerical results indeed show little change of the ratio as Q 2 is varied over a wide 
range of values for both singlet and non-singlet quarks and for gluons. 

We evaluated the influence of medium modification due to multiple scattering 
and induced gluon radiation on the above ratio in a cold nuclear medium. The results 
were found to be consistent, both with the experiment and also the vacuum evolution 
of the ratio. Calculations for the modification in a deconfined medium are in progress. 

References 

[1] M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 (1990); X. N. Wang and M. Gyulassy, 

Phys. Rev. Lett. 68, 1480 (1992). 
[2] C. Adlcr et al. [STAR Collaboration], Phys. Rev. Lett. 90, 082302 (2003). 

[3] J. C. Collins, D. E. Soper, and G. Sterman, in Perturbative Quantum Chromodynamics edited 

by A. Mueller, (World Scientific 1989), and references therein. 
[4] V. N. Gribov, and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); Yu. L. Dokshitzer, Sov. 

Phys. JETP 46, 641 (1977); G. Altarelli, and G. Parisi, Nucl. Phys. B126, 298 (1977). 
[5] A. Majumder, and X. N. Wang, Phys. Rev. D. 70, 014007 (2004). 
[6] A. H. Mueller, Phys. Rev. D. 18, 3705 (1978). 
[7] A. Majumder, and X. N. Wang, to appear, LBNL-51699 (2004). 
[8] P. Di Nczza for the HERMES collaboration, J. Phys. G30, S783 (2004). 

[9] X. F. Guo and X. N. Wang, Phys. Rev. Lett. 85, 3591 (2000); X. N. Wang and X. F. Guo, 
Nucl. Phys. A696, 788 (2001). 
[10] A. Majumder, and X. N. Wang, in preparation. 



